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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) 
UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to metric (SI) 


units as follows: 


Multiply by To obtain 
inches 25.4 millimeters 
2.54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144. meters 
square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.1745 radians 
Fahrenheit degrees 3/9 Celsius degrees or Kelvins' 


To obtain Celsius (C) temperature readings from Fahrenheit (F) readings, use formula: C = (5/9) (F — 32). 
To obtain Kelvin (K) readings, use forumla: K = (5/9) (F — 32) + 273.15. 


STORM SURGE SIMULATION IN TRANSFORMED COORDINATES 
Volume I. Theory and Application 


by 


John J. Wanstrath, Robert &. Whittaker, 
Robert O. Retd, and Andrew C. Vastano 


I. INTRODUCTION 
1. Background. 


Storm surges are transient fluctuations in the sea level induced 
by atmospheric disturbances, notably those due to extra-tropical 
storms and hurricanes and to a less frequent extent pressure jumps 
associated with line squalls. The rise of the water and circulation 
caused by a hurricane can be considerable and is of special practi- 
cal importance with respect to loss of lives and property not only 
adjacent to the coast but also well inland. Statistical studies of 
hurricanes of record provide a means of predicting the surge height 
along an open coast. The empirical formulas developed from these 
studies relate the maximum expected surge height to meteorological 
parameters and effective coastal bathymetry (Donn, 1958). However, 
all such studies do not provide the time history or even the time 
scale of the forcing function which is necessary as input for bay- 
response studies (Reid and Bodine, 1968). 


More recently, time-dependent models based upon the physics of 
the storm surge phenomena have been developed to study the generation 
or modification of the surge as it leaves deep water and moves over 
the Continental Slope and Shelf. These models, like the one proposed 
herein, involve the vertically integrated equations of motion and 
mass continuity. The greatest difficulty in utilizing these models 
has been the manner in which the shoreline has been portrayed and the 
application of realistic boundary conditions at the specified shore. 
Jelesnianski (1967, 1972) takes the shoreline as a vertical plane of 
infinite height, thereby facilitating the mathematical representa- 
tion of the shore boundary. More general portrayal of the shoreline 
is achieved by the schemes of Miyazaki (1963) and Platzman (1963) in 
which the coastline is represented as a series of straight-line seg- 
ments connected at right angles. Specifying the shoreline in this 
stairstep manner results in greater numerical programing complexity 
for the shore boundary condition in that the algorithm must possess 
the ability to search for and substantiate the location of land. A 
more serious objection is that this approximation may inject spurious 
oscillations into the calculations. This adverse effect is more than 
academic since the concentration of energy is near its maximum at the 
coast, and further, this is precisely the region from which water 
level observations are usually available for model verification. 


Clearly, these difficulties are avoided or minimized if a curvi- 
linear coordinate system is employed which will map the interior 
region bounded by the actual coastline, the seaward boundary (taken 
as the 180-meter depth contour) and two lateral straight lines. Thus, 
the transformation is such that the coast, as well as the 180-meter 
isobath, is mapped as a straight line in the image plane. In order 
to minimize the associated modification of the differential equations 
governing the storm surge, an orthogonal curvilinear system is 
employed in the present study, the desired mapping being carried out 
via a conformal transformation. 


2. Objectives. 


The objectives are twofold: (a) to develop a method which will 
conformally map the interior region bounded by two arbitrarily shaped 
curves each beginning and ending on parallel straight lines drawn 
between the two curves; (b) to develop and verify a numerical model 
of a hurricane-induced surge on the Continental Shelf employing the 
transformed equations of motion in the orthogonal curvilinear 
coordinate system. 


II. CONFORMAL MAPPING 
1. Development. 


The selection of appropriate conformal mapping and error functions 
is of principal concern in achieving the first stated objective. We 
desire to map a region R of the Z-plane into a rectangle in the c- 
plane (Figure 1), in which the seaward boundary and coastline curves 
are specifically transformed into the image plane as constant values 
of n . Furthermore, the curves in the negative x region are re- 
quired to be the mirror image of those in the positive x region. 

As an artifice to assure that the lateral boundaries of the mapped 
region represent straight parallel lines normal to the shoreline and 
the seaward boundary, symmetry about one of these boundaries (x = 0) 
is imposed and the whole range in x (-A to A) is considered to be one 
wavelength of a periodic function. Thus only the range 0 < x < i 

in the top panel of Figure 1 corresponds to the real region of the 
shelf. 


The conformal mapping relation is taken in the form: 
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An appropriate form of the transformation for the mapping considered 
in Figure 1 is: 
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where the coefficients A, , B, , C, , and D, are real quantities. 
However, the symmetry condition about x = 0 (or € = 0) requires that 
all A, and D, be zero. The resulting relations for x and y in 
terms of & and pn are: 


N 
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n=1 
and 
N 
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The condition that the range of x and € be the same requires that 
C, = 1. The remaining N values of C, and N+l values of B, are 
determined by matching the coastal and seaward boundary curves at 

n = + 6, respectively, 8 also being a parameter to be determined. 


The coordinates X,Y of the given coast or seaward boundary 
curves are specified parametrically in terms of arc length measured 
along each curve from some fixed point. The functions NSU AY ote. 
X° , and Y° are stngle-valued functions of this parameter where the 
superscript s or c represents the seaward boundary or coastline 
curve, respectively, This property is essential since a Fourier 
series-type representation is employed in determining the coefficients. 
The problem is to determine 8 , Bo , and the set of coefficients, 
By and Cy, , for a given N such that equations (9) and (10) give a 
best fit to the given curves, in the sense of minimizing an appropriate 
mean square error function, Since the specified curves are not known 
directly in terms of €& , but, rather in terms of arc length, the 
bicurve fitting equations require an iterative process starting from 
some initial estimate of arc length in terms of €& for each curve. 


It is found that a convergent iterative procedure results if one 
chooses the following as the i+l approximation of the coefficients: 
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where i=1,2... denotes the iteration number, A is the transform 
generated arc length and 
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The W terms are weighting factors used to distribute the error pro- 
portionally between the components of the transform generated and 
specified curves. 


It can be shown that the above relations result from minimizing 
the following error function: 
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The transform generated arc length and weighting factors are deter- 
mined relative to the results of the previous iteration by: 
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where E is the average value of equations (21) through (24). The 
arc length functions are given by: 
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Or OM ge < A and y° and yo are prorating factors taken such 
that A(A,n) is exactly the known length of the given curve. Two 
levels of approximation are involved in the error function. The ith 
approximation of these quantities is involved in the interpolation of 
the given curves. This approximation is mandatory in deriving a 
deterministic form for these quantities relative to the i+l level. 


The iterative procedure is initiated by setting all coefficients 
to zero, all weighting factors to unity, and equations (29) and (30) 
are given by yS&t and yet , respectively. The first 
estimates of 8 , BS and the Fourier-type coefficients 


are determined from equations (11) through (14), respectively. Hav- 
ing the first approximation for these quantities, the transform- 
generated arc length for each curve and weight factors may be deter- 
mined thus permitting a second approximation. This procedure is 
repeated for given N until 8 , B, and the transform coefficients 
have converged as indicated by successive values of the error func- 
tion. The error function approaches a constant value governed by N 
and the inherent errors in estimating the various integrals by 
numerical methods. The trapezoidal rule employing M values of € 
at equal increments over the interval 0 < € <A is used to estimate 
all integrals. Given the value of €& , the appropriate values of 

XS , YS , X° and Y° are computed by an interpolational subroutine 
called CURV, which was obtained from Dr. Alan Cline, National Center 
for Atmospheric Research. This routine is based on a numerical 
analogy to a spline under tension (for example, see Schweikert, 
1966). In all applications of the conformal mapping procedure, 

the total number of terms of the series, N , is taken as M or 
smaller. The limitation on N is imposed since, for a determinate 
system, the total number of Fourier-type coefficients cannot exceed 
the number of data points. The selection of N is governed by the 
desanedudegreeyot tation ches functionss Xen) Yoo se Xo anda Your 
The procedure outlined above follows the rudiments of a conformal 
mapping technique developed by Reid and Vastano (1966). 


2. Applications. 


Three regions of the Continental Shelf of the Gulf of Mexico and 
two regions of the eastern seaboard of the United States have been 
mapped. These regions are: 


(a) western gulf coast: Laguna Madre, Mexico, to Marsh 
Island, Louisiana; 

(b) central gulf coast: Matagorda Bay, Texas, to Timbalier 
Bay, Louisiana; 

(c) eastern gulf coast: Atchafalaya Bay, Louisiana, to 
Apalachee Bay, Florida; 


(d) lower east coast: Cape Kennedy, Florida, to Pamlico 
Sound, North Carolina; 
(e) upper east coast: Pamlico Sound, North Carolina, to 


Penobscot Bay, Maine; 


and are shown in Figures 2 through 6. The top and bottom curve in 
each figure represents the coastline and the 180-meter isobath, 
respectively. The latter is taken to correspond approximately to 
the shelf break, seaward of which the depths increase abruptly; the 
180-meter contour may not always be the most appropriate one for the 
outer limit of the Continental Shelf. The chart scale, orientation 
of the region, and location of cities and bays along the coast are 
indicated on each figure. The x's represent the location where 
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the coast and the 180-meter contour were digitized from an overlay 
of bathymetric charts for the regions concerned. The smooth curve is 
the result of applying the numerical spline interpolation routine to 
these points. In reading the discrete positions, both the coastline 
and the 180-meter isobath have been smoothed subjectively to suppress 
variations with scale lengths less than the grid scale ultimately 
employed in the surge calculations. For example, narrow entrances 

to bays are replaced by a fictitious coastline across the entrance, 
and the cusp-shaped features are smoothed. 


A method was adopted which minimizes the amount of computer time 
required to solve equations (9) and (10) by hastening the convergence 
of the iterative procedure. The method consisted of increasing N 
at selected iteration intervals. The iterative procedure used to 
determine the transformation coefficients was terminated if the con- 
vergence criterion was satisfied or if the available computer storage 
was exceeded because of increasing N , The convergence criterion 
was a mean variance of less than 1 square kilometer between the 
transform-generated curves and that specified, However, for develop- 
ment purposes, the iterative procedure for the western and central 
gulf coast and lower east coast areas was continued beyond this 
criterion to obtain a better fit, 


It is convenient to discuss the mapping of the western and central 
gulf coast regions at the end of this section since additional testing 
of the solution to the mapping equations was performed with these 
areas. Figure 4 shows the shelf region for the eastern gulf coast 
where the Mississippi Delta is shown as the shaded area. Considera- 
tion of the numerical time step for the surge algorithm was the prin- 
ciple reason why the coastline to be mapped did not follow the delta. 
The reduction in the time step if the actual coastline had been 
followed would, probably, have been at least tenfold with respect to 
the one used in the simulation of the storm surge induced by Hurri- 
cane Camille (Section V). Since the delta, or really the levee 
adjacent to the Mississippi River, has a controlling influence on 
the circulation and surge caused by the hurricane, this geographical 
feature was included in the surge model as a wall protruding from the 
coast. The fit of the transform-generated coast and seaward boundary — 
curves with respect to that specified after 1, 20, 40, 80, and 160 
iterations is shown in Figures 7 through 11, respectively. The suc- 
cessively better agreement of the mapped curves with respect to that 
specified is obvious from the convergence table for the gulf coast 
region of Atchafalaya Bay to Apalachee Bay (Table 1). The table 
shows the variance and convergence behavior of selected transforma- 
tion coefficients as the iteration and N increases. The coast 
variance in X , the third column of Table 1, is calculated by: 


M : 
var PS Gio) ae. P (31) 
jel 
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Similar expressions are employed for the other variances. Appendix 
A contains the transformation coefficients used to produce Figure 11, 


The successful application of the mapping equations to the other 
shelf regions was accomplished. Tables 2 through 5 clearly indicate 
the successive convergence of the mapped curves with respect to that 
specified for the lower and upper east coast and the western and 
central gulf coast, respectively. Figures 12 through 15 show the 
fit of the transform-generated curves and that specified after the 
last iteration for the above regions, Appendixes B and C contain 
the mapping coefficients used to produce Figures 12 and 13, respec- 
tively. The bottom part of Tables 4 and 5 show the results of addi- 
tional testing of the conformal mapping equations (9) and (10) with 
the western and central gulf coast regions. Another, less general, 
solution of equations (9) and (10) is possible if one minimizes the 
least square error function defined only in terms of the Y inte- 
grals. This solution for the transformation coefficients, hereafter 
referred to as an alternate solution, may be obtained from equations 
(13) through (19) with we = We = 0. The testing procedure was to 
continue the iterative process as outlined in the previous section 
with the initial approximation for the coefficients being those 
values determined from the 80th iteration. The alternate solution 
as applied to the western and central gulf coast regions was stable 
and, moreover, provided a better fit with N = 110 than the more 
general one. The fit of the mapped curves with respect to that 
specified after 1 and 40 (or, 81 and 120) iterations with the alter- 
nate solution is shown in Figures 16 and 17, respectively, for the 
western gulf coast and Figures 18 and 19, respectively, for the cen- 
tral gulf coast. Appendixes D and E contain the transformation co- 
efficients used to produce Figures 17 and 19, respectively. In test- 
ing with the other three shelf regions, the alternate solution was 
nonconvergent in that the successive values of the error function 
do not decrease or approach a constant as outlined in the previous 
section. The reason for this behavior has not been investigated, 


The total variance of x and y froma linear transform of the 
curvilinear coordinates may be defined as: 


Ge = ae + oh (32) 
where 
BS A 
ea laa 
eh) TDN iS BJP cle alg 3 (33) 
-B Oo 
and 
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Bin 
ges aa | | ly - + n)]2 dé dn . (34) 
Bano 


These expressions are measures of the “"curvilinearity" of the trans- 
forms x(&,n) and v(é&,n) . From equations (9) and (10), it can be 
shown that 


N 
or = ) Cis 5 (35) 
n=1 
where the contribution to the total curvilinearity variance due to 
the nt” term is: 


22 2 Q 
(BL + Co ) sinh 2nkg 


Dales 
ons 4nkg F Ce) 


The value of oe ee is a more meaningful measure of the relative 
A Gc of the various terms in the transformation than the raw 
coefficients. The curvilinearity variance spectrum given 

on 22 as a function of n_ for the final transform of the five 
ae is shown in Figure 20. The values of o% and o/d for’ the 
five regions are indicated in the figure. Each of the curves indi- 
cates a general inverse power law trend of about 4th degree. The 
variability of the spectrum in all but the middle curve is more pro- 
nounced in the higher harmonics than in the lower. 


3. Limitations. 


Testing of the conformal mapping equations was conducted with 
coastlines that varied in the degree of smoothness and regions that 
varied in i . Another version of the coastline from Cape Kennedy 
to Pamlico Sound is shown in Figure 21. Comparing Figure 5 with this 
version, it can be seen that there are two small differences in the 
form of the coastline. Figure 21 shows the narrowing of the Con- 
tinental Shelf in the area south of Cape Kennedy. Another differ- 
ence is the slightly sharper vertex at Cape Kennedy, Cape Fear, and 

Cape Lookout. The (general) solution of the mapping equations proved 
unstable for this situation. It is not known which change in the 
coastline permitted the convergence of the iterative procedure for 
Figure 5. 


One test of the mapping equations was conducted on the region 


from Matagorda Bay to Apalachee Bay (Figure 22). This region is the 
result of joining Figures 3 and 4 and represents in this study, the 
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Figure 20. The “curvilinearity" variance spectra for the final 
transform of the five mapped regions. 
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longest horizontal extent of a region to be mapped. The mean vari- 
ance between the mapped curves and that specified after 120 itera- 
tions with n=140 was 4.4 square kilometers, The level of vari- 
ance achieved indicates that the bicurve fitting equations are 
marginally successful in mapping this enlarged region, The change 
in the variance level over the last 10 iterations and limits on 
computer time were the bases for terminating the iterative procedure. 
The convergence criteria can be satisfied either by a shorter hori- 
zontal extent or smoother curves for a longer region. 


49 


III. THE STORM SURGE EQUATIONS IN THE 
SHELF COORDINATE SYSTEM 


1. Stretched Shelf Coordinate System. 


Consider the transform-generated coastline and seaward boundary 
curves as shown in Figure 23. The orthogonal curvilinear mesh asso- 
ciated with the coordinates (&,n) is designated the shelf coordt- 
nate system. Although the surge calculations can be performed in 
the shelf coordinate system (Figure 23b), a problem remains in 
obtaining the desired spatial resolution with the fewest possible 
computational points. The computing grid employed (Figure 23c) is 
the result of a second transformation which preserves the orthogo- 
nality property but allows independent stretching of &€ andof n. 
The grid resulting from the second transformation is termed the 
stretched shelf coordinate system (S*,T*). 


The stretched shelf coordinate system is generated by indepen- 
dently transforming the & and n axes in the following manner: 


(a) Given the hurricane storm track, the location where 
the storm crosses the coast and the coastline con- 
figuration, the values of &€ along the coastline 
are determined which will produce a constant rela- 
tively fine increment of coastline arc length, Sp 
In this area of prime interest, line BC in Figure 
23(a), the constant increment of arc length is 
equal to AS* . However, regions AB and CD show 
that for the same AS* as above, there is a rela- 
tive expansion of the increment of the coastline 
arc length. The functional relationship between € 
aimcl $9 wise 


S* = 5S (ey) ? (37) 


where the expansion of S, with respect to S* 
is specified by an (arbitrary) analytical ex- 
pression of the form = A + B(S*)C where A , 
B, and C are constants. 


(b) For the shelf bathymetry along a particular centrally 
located isoline of & , the values of n are de- 
termined which will yield a constant change in the 
time, AT* , required for a long wave traveling at 
the local free wave celerity to proceed from seaward 
boundary to the coast. The long wave travel time, 

T , is calculated by: 
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(a) 
Re 
ves & 1 
| inns CTE 
| FERRE ES 
SO aee e, 
depth contour 
0 
0 errs ce ; 
(b) 
(c) 


Figure 23. Orthogonal curvilinear grid system in the Z-plane 


(a), in the ¢-plane (b), and in the computing 
grid (c). 


5! 


; (38) 
so) ) 


where) S57) 1s) chendisitancer a longmthestsolinevotame uD ytsmene 
local depth relative to mean sea level, and g is the acceleration 
due tg gravity. The relation between n and T* is given by: 


T* = TG ())) , (39) 


The incremental values of T* are determined from equation (38) 
subject to the (arbitrary) expansion relationship of T*(T). This 
relationship is a convenience which permits an additional degree of 
freedom in adjusting the relative spacing between isolines of n . 
However, if T* = T , the value of AT* is that which divides the 
total long wave travel time by an integer number of lines of nn. 

The selection of the final AT* is based upon a compromise for pro- 
viding adequate resolution for the wind field and the resulting surge 
with a minimum number of points. The form chosen for the expansion 
13, IF eo Ae B(T#)© where A, B, and C are constants. 


The stretched shelf coordinate system provides a grid system 
with a finer resolution near the coast than at the 180-meter iso- 
bath. The expansion curve, S*(Sp) , stretches the horizontal reach 
of the grid while maintaining a finer grid in the area of landfall 
of the hurricane. In this manner an economy is achieved in terms of 
the number of grid points required by the surge program, However, 
because the preferred expansion curves dictate the locations (in 
x,y space) of the depths required for the surge calculations, the 
depth field must be redefined for different pairs of stretching 
functions. Application of this procedure has been accomplished for 
the simulation of the storm surges caused by Hurricanes Carla, 
Camille, and Gracie. Figure 24 shows the transform-generated coast- 
line and seaward boundary curves for the western gulf coast area. 
Also shown is the shelf coordinate system (&€,n) which was employed 
in the Hurricane Carla surge simulation. The track of the hurricane 
is the dashline. The grid system extends over approximately 750 
kilometers of coastline, The stretched shelf coordinate system 
(Figure 25) is attained through the functional relationships speci- 
fied in Figure 26 for the S* axis and Figure 27 for the T* axis. 

The shelf coordinate system, the stretched shelf coordinate system, 
and similar transformation relationships are shown in Figures 28 
through 31, for the grid used to simulate the storm surge caused by 
Hurricane Camille, and in Figures 32 through 35 for the grid asso- 
ciated with the Hurricane Gracie storm surge. 
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Shelf coordinate system (&£,n) for Hurricane Carla surge simulation. 
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Figure 29. 
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Shelf coordinate system (£,n) for Hurricane Gracie surge simulation. 


Figure 32. 
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2. Storm Surge Equations. 


The vertically integrated form of the quasi-linear long wave 
equations in a Cartesian system are well known (e.g., Welander, 
1961). The appropriate forms for storm surge computation in the 
stretched shelf coordinate system are as follows: 


dQ 
HORS Ti) Br 2D OR as iN i 
Aite fy Qe ts Fu aS* (H Hp) i Tox Tox > (40) 
dQ 
T* gD 3 i . 
yr See ey pre BS) = Gs 9 ore ~ (en 
and 
SHE ual ates 13 \ 
ara ap a as (FQcx) tp Ser Glee)! = (0), (42) 


where Q is the volume transport per unit width, t is the wind 
stress divided by the water density, o is the bottom resistance 
stress divided by water density, f is the Coriolis parameter, D 
is depth of water, H is the sea surface elevation relative to mean 
sea level, and Hp is the hydrostatic elevation of the sea surface 
corresponding to the atmospheric pressure anomaly. The term F “is 
a variable scale factor associated with the conformal part of the 
orthogonal curvilinear coordinate system. Specifically, 


F= (()2+ @7, (43) 


where x and y are given by the transformation equations. The 
terms wu and vy , are additional scale factors that relate to the 
transformation of &€ to S* and n to T* , respectively, and 
are given by: 


3E aS 
Me ee ae ve (44) 
ie 
and 
3S 
Gy 
y = -t goon eb M (45) 
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Furthermore, F*yv corresponds to the Jacobian of the transforma- 
tion in the sense that 
Nea fF Fi dss cite; (46) 
Re 
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where A is the area of a closed region R in the x,y plane 
whose mapped image is R* . 


The kinematic wind-stress components T., and 
lated to their x,y component counterparts Cr 5 
point by: 


Tre are re- 
Fo) at a given 


uae = cos 6 + fe sin 6 , (47) 
and 
Tp. = “Tt, Sin @ + Sy cos 6 , (48) 
where 
d = (Oy/OE 
Om=mtan Gama) (49) 


(50) 


where p, is the air density, pw is the water density, Cy is 
nondimensional drag coefficient, and W,9 is the windspeed at an 
elevation of 10 meters above the water surface. The value of K 

(Reid and Bodine, 1968) is taken as; 


Ky » see Wag, < oO moses 
Ke -l 
Ky, + [1.0 - 7.0/Wy9]? Ko , if Wyo > 7.0 m-s 
MMe Ky ds lei sa0°S and kK as 2.5 % 1076, 


The values of K, and Ky used by Whitaker, Reid, and Vastano 
(1973) for Lake Okeechobee were found to be too large for the shelf 
surge computations investigated in this study. 


The form of the bed resistance terms are: 


K,Q 
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and 


Ke 


T* = D2 Qn ’ (53) 


0, 
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where 


2 2435 
OQ lle, 7 Grell” 6 (54) 
and Ky is a nondimensional drag coefficient taken as 2.5 x 1073. 


3. Numerical Algorithm. 


The numerical analogs of equations (40), (41), and (42) are based 
upon centered difference approximations of all terms (see App. F). 
The algorithm treats the time dependency explicitly and employs a 
computing lattice as shown in Figure 36 in which the transports, Q., 
and Qrsx are computed at the same location but are staggered in S 
time and space with respect to tne water level anomaly. This scheme 
facilitates a simple representation of the Coriolis and bottom stress 
terms in the difference equations. The surge model allows for vari- 
able bathymetry and arbitrary coastlines which are represented as a 
series of connecting straight line segments situated along lines of 
constant S* or T* . For the simulation of the Hurricane Camille 
storm surge, the arbitrary coastline feature of the algorithm is 
mandatory to delineate the delta. However, in the more usual appli- 
cation, the coast is a straight line in the computing grid. 


The difference equations for Qcx, Qrx, and H at interior 
points of the computing grid are given by: 


i 


OL Cipal) os (Gi Go Ne Gal (Gr 2 G2 Maye) } (55) 
S 


Qre(isdomtl) = (G; G3 - f at GiiG 2 & hey) » (56) 


and 


H(i,j,nt+l) = H(i,j,n-1) 


f At IP((aleri 59) Qo. Gt1,j,n) ~ F-1,j) Qo. G-1,j,n) 
BGI) u(a) AS* 
Es Giewital) Qpa G2, j+1,n) a (ab, 9} JL) Qpx G,j-1,n) 

i 0G) aT oe 


(ines als ' 5) 5 Gicl m iinchices Gymesss wag SP 5 WY 5 eine! teams 
coordinates, respectively, and At is the numerical time step. The 
quantities G, , G» , and G3 are given by: 


G) = 142K, at QGi,j,n-1)/D7 , (58) 
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Figure 36. Scheme for computed variables. 
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u(i) ss) AS* lil (Gaisell) 5in)) H,(i+1,j,n) 
- H(i-l,j,n) + HyG-1,3.n)) + 2 At Tox ; (59) 

G3 = Qrx Gs j.n-1) ae aerAN Qo, 4,j,n-1) 

4 g At D Ate Zi Gah 
vj) FA Ge) AT* HG, 51.0) Hp (4, j+1,n) 
- H(i,j-l,n) + Hy@i.j-1.0)| i 2D IME Tipe > (60) 

where 
D = A D(i+1,j,n) + Y(t Ik 5 3 50) at D) (Gis aerdl 5101) ty DG,j-1.)| * (61) 


In order to maintain numerical stability, the time step must 
satisfy the condition: 


Fuv AS* AT*) 5 
At <[ Ane ye (62) 


where 


s 


G, = [(@m AS3)2 & Gy areyzy" , (63) 


and the minimum value of the right-hand side of relation (26) is 
implied. Thus, a search of the grid is required for a proper selection 
of At . However, the conditions along the seaward boundary usually 
limit the time step because of the deeper water. 


The values of »p , v=, F , and 6 given by equations (44), (45), 
(46) and (49), are required at the appropriate computational grid 
points. A numerical spline under tension is employed to interpolate 
& given the IM wallues oF S (S*) . Wne scale Accor, mM 5 1S 
determined by: P 


ua) = BGs) eb (64) 
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where €(0) and €(IM+1) are ascertained by expanding the computing 
grid one increment of AS* to the left and right, respectively. In 
a similar fashion, the JM values of 1 may be determined and the 
scale factor v given by: 


vie) = MUD = alien) 


2 (ST 2 ed) 
J =& 2, e806 JME. 
Since the computing grid may not be extended landward or toward deeper 
water, the approximation at j = 1 and j = JM is made that: 
wa) = In@) = a@ylfare (66) 
and 
v(JM) = [n(JM) - n(JM-1)]/AT* . (67) 


The area in the x,y plane of each grid block, including those 
from expanding the grid in the & direction, is determined by numer- 
ical procedure from the following relation: 


Me ei s 
Lote 2 am | 
A(@i,j) = | 2D lea. (68) 
0g OE 
Tee oS 
J 1-1 
i 5 1,2 IM+1 , 
1,2 1,2.¢ JM 


The scale factor associated with the conformal mapping is determined 
from equation (46) at intersectional points, excluding the coast and 
seaward boundary, by the expression: 


z I. 
Be) < Mas) ee A(t jel) = AGieil 7) & Atl, je) | - 
>J AINSI icant GN) aol Gin) 


(69) 
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and for grid points along j = 1 or JM by the approximation: 
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FL.3) = la-ag* at Gp) wG) : (70) 


th eZ scteh MLM ied 


The angle, 6 , relating the orientation of the stretched shelf 
coordinate system to the x,y plane is determined at grid points, 
excluding j = 1 and JM, by the following smoothing procedure: 


ated, | a ft rita é : oa Tae 
(53) - =|6(4-1,3) + dG gen) = d@an.aye d(i,j+1) + 48.3] 


(71) 
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j= 2,3 °° JM-1, 


where e! is given by equation (44). Along the sea and coastline 
boundaries we have 


e(i,1) = |e (i-1,1) + 20'(4,2) + g! Gsugi) s 40" Gi,1)|  &® 
4 21,2 ose ing | 
and 
aC ,an = de" Gta 3s eA NS 8 Galote Asta 
f= 1g) 392) iM 7%) 


4. Boundary Conditions. 


A wall condition is employed in the surge simulation along the 
coast while the surface elevation anomaly is placed in equilibrium 
with the atmospheric pressure along the seaward boundary. Vanishing 
normal derivatives of transport are specified on the lateral open 
boundaries. This condition is used by Jelesnianski (1965, 1966) and 
Forristall (1974), although it may not be the most desirable (Reid, 
O75) < 


In the x,y plane the coast is curved making the wall (coast) 
and lateral boundary conditions difficult to apply in a rectilinear 
erid system. However, the stretched shelf coordinate system repre- 
sents the coast as: 


Qe sO. (74) 


or, the analog is simply, 


Ta 


Qra(i,JM,n+1) = 0, (75) 
1 S158 ++- TM, 

Thus, total reflection is assured at the boundary. The flux, Qgox, 
along the coastline is calculated from equation (55) with D taken 
as the two point average of the local fluid depths along the coast. 
The local depths ranged from 1.0 to 2.3 meters, depending on the sea- 
ward bottom slope and surrounding elevations. The water elevations 
along the coast are computed from the continuity relation as given 
by equation (57) with the following substitution: 


F(i,JM+1) Qo, (i,JMt1,n-1) = -F(i,JM-1) Q.,(4,JM-1,n-1) , (76) 


1S Aya ooo EZ 
This is an artifact consistent with total reflection. 

For the simulation of the Hurricane Camille storm surge, the 
normal routine of the surge program was interrupted at those grid 
points representing the protruding Mississippi Delta. Along this 
part of the coast (the solid heavy line in Figures 28 and 29), the 
normal flux was set to zero and the tangential flux was determined 
from either equation (55) or equation (56) with the Coriolis term 
vanishing. The continuity relation is altered depending on the 
orientation of the boundary to be consistent with total reflection. 

The open deep sea boundary condition is: 

H@, l5n+l))) = H,G,1,n+1) i (7) 
io 24 coo MME 
and 
H(i,2,n+l) = [HG@=1,1,n+1) * H@i+l,1,n+1) 
+ H(i-1,3,n+1) + H(i+1,3,n+1) ]/4 , (78) 
1 8 O95 999 NEP 
Specifying the seaward boundary in this sawtooth form obviates the 
calculations of Qcx and Qrx along the boundary since they are not 


required for computations at interior points. 


The lateral open boundary condition requires the normal gradient 
of the S*-directed transport to vanish, 


2 


9Qox 
Sem Tee (79) 


This condition implies that along the left side of the grid, 
Qox 1,5 ,ntl) = Qo (S,j,nt1) , (80) 


==) Ove 
and along the right side, 
Qox(IM,j.m+1) = Q.,(IM-2,j,n+1) , (81) 
Te Sy OO HM ae 


The S*-directed transports for even j and i=2 or IM-1 are deter- 
mined by the average of the two neighboring interior values. Addi- 
tionally, the T*-directed transports are required along i=2 or IM-1 
and are calculated from equation (56) where D is again the average 
of the local fluid depth along the boundary, 


Special computations are required at seaward corner points (2,2) 
and (3,1). Simultaneous equations are solved for Qpx(2,2,n+1) and 
Qox(3,1,n+1) with the approximation that Qypx(3,1,n+l1) is the aver- 
age of the two neighboring interior values. Similar expressions are 
used at the right-hand seaward corner points. 


Other conditions at the seaward and lateral boundaries were 
experimented with for the Hurricane Carla simulations. These in- 
cluded radiational conditions of the type discussed by Reid and 
Bodine (1968) and Reid (1975). The main differences, as anticipated, 
were close to the lateral boundaries where the "flow through" condi- 
tion in equation (79) gave a more realistic rendition of the long- 
shore flow. 


The basic numerical model was used successfully by Alvarez 
(1973). Further testing of the algorithm, where in some cases an 
analytical solution was possible for comparison purposes, is pre- 
sented in Appendix G. 


5. Wind and Pressure Fields. 


Of major importance in simulating the storm surge is the accurate 
portrayal of the hurricane wind field on the computing grid. More- 
over, the representation of the wind must be time-dependent. 


The Hydrometeorological Section of the National Weather Service 
(NWS) provided charts of the surface winds at 10 meters above mean 


>) 


sea level and of the barometric pressure for Hurricanes Carla and 
Camille. One method for representing the wind field is to digitize 
the above charts for the hurricane in question. Usually the wind and 
pressure data depicting the storm are sampled in time at particular 
grid points and interpolated to provide the necessary input to the 
surge model. However, such input data are tedious and laborious to 
obtain and do not guarantee better quality of the input than can be 
obtained by analytical representations of the forcing fields. It is 
this latter approach which is used in this study. 


For a given hurricane, the following parameters are sampled in 
time from the NWS charts: position, forward speed, central pressure, 
radius vector (relative to the storm movement) to maximum winds, and 
the maximum winds. With an analytical representation of the surface 
wind and pressure fields the above parameters need only to be inter- 
polated in time. 


The wind field representation prior to consideration of the land 
influence as given by Jelesnianski (1965) is employed in this study. 
The x,y-wind components for a stationary storm are: 


V 
apoksiie. 27 iat 
Woe = Tp [ (x x) sin 9 (y Yn cos ] 
and F(r,) > (82) 
VR 
Uy = z, | + (x-x, ) cos > - (y-y,) sin ¢ | 
where 
2 45 i 
3 2 (eee) © Goa de 1° (83) 
Vp is the maximum wind, 9» is the ingress angle, and F(r,) is 
BYO I/2 « 
(r,,/R,) ge ig. S RB or (R/T) if ry > R- Distance from 


the storm center Oy,» Yn to the region of maximum winds is RX é 


The translation of the storm provides an alteration in the wind field 
which is carried out following the method of Jelesnianski (1965). 
This involves the vector addition to the above field of a supplemen- 
tal velocity whose direction is parallel to that of the storm and 
whose magnitude depends only upon xy, . The resulting modified 
values of V_ andV_ are used to compute the wind-stress 
components: y 
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This constitutes the symmetric wind-stress field. ‘The stress compo- 
nents in the stretched shelf coordinate system are determined by 
applying equation (84) in equations (47) and (48). 


The surface atmospheric pressure field associated with the hurri- 
cane is given by: 
“Ra! Th 


Peps @ Pye ; (85) 


where P, is the central pressure and P, is the far field pressure. 
The term Hp in equations (40) and (41) 1s equivalent to 
(= EUR Se 


The symmetric surface stress field does not reflect the influ- 
ence of land. Without modeling this influence the analytically 
determined wind field near the coast would not be consistent with that 
from the NWS. A systematic procedure (App. H) is employed to alter 
the symmetric wind field such that it portrays the effect of land. 
These winds are referred to as deformed. Mr. Thomas Reid, Texas A&M 
University, Department of Oceanography, is responsible for the basic 
development of this wind model (unpublished manuscript). The pressure 
field was not altered in the nearshore region to conform with the 
wind field, 


Figure 37(a) represents the symmetric Hurricane Carla wind field 
at 1600 Greenwich mean time (G.m.t.), 11 September 1961 (approxi- 
mately 4 hours before the storm crossed the coast), Isovels are 
shown in meters per second. The analytically deformed wind field at 
this time is shown in Figure 37(b). The influence of land on the wind 
field is evidence even at a time when the storm is approximately 55 
kilometers from the coast. Figures 38 and 39 show the symmetric (a) 
and deformed (b) winds at 2000 G.m.t.,11 September, and 0000 G.m.t., 
12 September, respectively, The rotation and reduction of the wind 
vector in the nearshore region to reflect the land influence is 
illustrated in comparing Figure 40(a) for the symmetric case with 
Figure 40(b) for the deformed wind. These figures correspond to the 
contours shown in Figures 38(a) and 38(b), respectively, The arrows 
represent the wind vector placed such that the tail is at the com- 
putational point. The pressure fields, as determined from equation 
(85), for Hurricane Carla at the above times are shown in Figure 
41(a,b,c) where the contours are in millibars, Although the pattern 
of isovels reflecting the land influence is in good agreement with 
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Time = 46.0 (hrs.) 
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Figure 37. Hurricane Carla symmetric (a) and deformed (b) 
wind fields at 1600 G.m.t., 11 September 1961; 
isovels in meters per second. 
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Time = 50.0 (hrs.) 


Figure 38. Hurricane Carla symmetric (a) and deformed (b) 
wind fields at 2000 G.m.t., 11 September 1961; 
isovels in meters per second. 
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Time = 54.0 (hrs.) 


Figure 39. Hurricane Carla symmetric (a) and deformed (b) 
wind fields at 0000 G.m.t., 12 September 1961; 
isovels in meters per second. 
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Time = 50.0 (hrs.) 


(a) 


(b) 


Figure 40, Hurricane Carla winds at selected points at 
2000 G.m.t., 11 September 1961 for the 
symmetric (a) and deformed (b) representations. 
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Time = 46.0 (hrs.) 


(b) 


9 Time = 54.0 (hrs.) 


(c) 


Figure 41. Hurricane Carla atmospheric pressure on 11 Sep- 
tember 1961 at 1600 (a), at 2900 (b), and 0000 
G.m.t., 12 September (c); contours in millibars. 


80 


those provided by NWS, the problem of specifying the proper deforma- 
tion is not completely resolved. 


The observed high winds which remain along the coast after the 
storm has proceeded inland are attained by setting R, to be the 
distance the storm center is from that coastal area of persistent 
high winds. This procedure was also followed in modeling the winds 
from Hurricane Camille. However, R, in the pressure expression 
was not allowed to increase in this manner. 


The hydrographs at selected grid points along the coast are 
presented in Figure 42 illustrating the improved agreement with 
respect to the observed water levels in using a deformed wind 
forcing function rather than a symmetric one, The observed water 
levels, corrected for the astronomical tide, are indicated by 
squares. The computed water levels have been raised for comparison 
purposes to correspond to the local sea level at the start of the 
computations (1800 G.m.t., 9 September). The solid line represents 
the computed surge employing the deformed wind and the dashline 
shows the results with a symmetric wind. The surge computations 
were formed using identical boundary conditions. 


The symmetric (a) and deformed (b) winds for Hurricane Camille at 
0000 G.m.t. 18 August 1969 (approximately 4 hours before the storm 
crossed the coast), 0400 G.m.t. and 0800 G.m.t. on the same date 
are shown in Figures 43 through 45. The rotation and reduction of 
the wind stress for this storm are different from that of Hurricane 
Carla because of the differences in the coastal configuration. 
Figure 46 (a,b) shows the wind vectors for the symmetric and deformed 
representation of Hurricane Camille at the time the storm proceeded 
inland. The pressure fields at the above times are presented in 
Figure 47(a,b,c). The wind fields associated with Hurricane Gracie 
were not deformed because the documentation was not amenable to the 
deformation procedure. Also, Rh was not increased in the manner 
specified above after the storm crossed the coast. The symmetric 
winds at 1200 G.m.t., 29 September 1959, 1600 G.m.t. (the time the 
storm proceeded inland), and 2000 G.m.t. on the same date are shown 
in Figure 48, Figure 49 presents the pressure fields at the above 
times. 
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Observed (squares) and computed water levels 
using a deformed wind (solid) and a symmetric 
wind (dashed) for Hurricane Carla. 
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Time = 24.0 (hrs.) 
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Figure 43. Hurricane Camille symmetric (a) and deformed (b) 


wind fields at 0000 G.m.t., 18 August 1969; 
isovels in meters per second. 
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Time = 28.0 (hrs.) 


Figure 44. Hurricane Camille symmetric (a) and deformed 
(b) wind fields at 0400 G.m.t., 18 August 
1969; isovels in meters per second. 
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(b) 


Figure 45. Hurricane Camille symmetric (a) and deformed (b) wind 
fields at 0800 G.m.t., 18 August 1969; isovels in 
meters per second. 
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Figure 46. Hurricane Camille winds at selected points at 0400 
G.m.t., 18 August 1969 for the symmetric (a) and 
deformed (b) representations. 
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Time = 24.0 (hrs.) 
(a) 


| Time = 28.0 (hrs.) 


(b) 


(c) 


Figure 47. Hurricane Camille atmospheric pressure on 18 August 
1969 at 0000 (a), at 0400 (b), and at 0800 G.m.t. 
(c); contours in millibars. : 
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Time = 24.0 (hrs.) 


Figure 48. Hurricane Gracie symmetric winds on 29 September 
1959 at 1200 (a), at 1600 (b), and at 2000 G.m.t. 
(c); isovels in meters per second. 
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Time = 24.0 (hrs.) 


(a) 


(b) 


9 Time = 32.0 (hrs.) 


(c) 


Figure 49. Hurricane Gracie atmospheric pressure oi 29 September 
1959 at 1200 (a), at 1600 (b), and at 2000 G.m.t. 
(c); contours in millibars. 
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IV. SIMULATION OF THE FREE WAVE IN AN ANNULUS 


1. Problem Statement. 


If the wind stress, bottom friction, atmospheric pressure, and 
the rotational effect of the earth are neglected, then the equations 
in Section III admit a simple analytic solution for free gravity 
waves of long wavelength. A study using curved boundaries will 
demonstrate the superiority of modeling the long gravity wave in 
orthogonal curvilinear coordinates over those models which employ 
rectilinear coordinates. Basically, the study involves comparing 
the numerical solution of a free-standing long wave in a 90° section 
of an annulus with an inner radius r, of 393 kilometers and an 
outer radius ry, equal to 2r, (Figure 50). The annulus is bounded 
on all sides by a vertical wall. The depth of the basin relative 
to the mean water level is assumed everywhere to be 40 meters. 


The analytical solution for the free-wave oscillations in a 
section of an annulus may be obtained with some modifications from 
Lamb (art. 191, 1932). The boundary conditions require that: 


oH L Oh he ap seh) Ele! see ae) (86) 
or 

and 
a7 = 0, at @ = 0° and @ = 0/2 . (87) 


The analytical solution for the given initial conditions is: 


H(r,6,t) = yi y A J (ge 32) = Talkan?}) ie) | 
ehh oe m,n n> m,n Y'(k r) 2% m,n 
cos né cos Gamat ‘ (88) 
m= 0,1 3 
n = 0,2 
where 
omn i Spi ved, (2) 


and for given n , kh 2 is the mth root of 
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Fisure 50. The annulus in polar coordinates (rT,€). 
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Ji(kr2) Y¥ (kr) - J) (kry) Yi(kra) = 0. (90) 


The terms J and Y are Bessel functions of the first and second 
kind, respectively. The order is given by the integer subscript and 
the superscript refers to differentiation with respect to the argu- 
ment. The term A in equation (88) is an arbitrary constant repre- 
senting the initial amplitude. The first azimuthal mode is for 

n= 2 and the lowest value of k that satisfies equation (90) is 
approximately 1.340/r,; meters-!. The period of oscillation for this 
mode (m = 0 , n = 2) is 25.85 hours. 


The numerical solution of the free wave in the quarter annulus 
is sought by performing the integration in two different computing 
grids. In one case the grid is rectilinear (Figure 51), and in the 
other a polar system is used. In the rectilinear system the outer 
and inner radii of the annulus (the light line in Figure 51) is simu- 
lated in a stairstep fashion. Considering a limitation on computer 
time and storage, an acceptable rendition of the curved boundaries 
is present in the Cartesian grid. Proper representation of the 
quarter-annulus requires a transport point at corners on the outer 
and tnner boundartes. Consequently, the rectilinear boundary is not 
symmetric about 1/4. The locations of nine hydrograph positions are 
indicated by small boxes. Although the computing grid is 43 by 43, 
only 1,052 points are used to represent the annulus. In this grid 
system, AS* and AT* are just Ax and Ay , respectively, with 
Ax = Ay = 19.65 kilometers. The maximum allowable time step as 
determined from equation (62) is approximately 700 seconds where F = 
u = v=1. Finally, the analog form of the long wave equations is 
obtained by setting those terms in equations (55) through (57) to 
zero which are neglected, and setting F , uy, and v to one. The 
transports, Qox 5 Qr* » are now aliases for Q. and Qe 5 
respectively. 


The numerical algorithm for the rectilinear grid system must be 
capable of determining if a computational point is interior, 
exterior, or on the boundary of the annulus. For exterior points, no 
calculations are performed. Grid points on the stairstep boundary 
require special attention subject to the condition of a wall. 
Furthermore, the numerical program must identify and apply the appro- 
priate wall condition depending on whether the point in question is 
located at a boundary corner (an H computation) or on a segment of 
the boundary. Clearly, extensive programing and computer time is 
required to accomplish this task. 


The other computing mesh which is used for obtaining the numerical 
solution of the free wave in the annulus is the polar (or stretched 
shelf coordinate) grid system (Figure 52). The transformation 
to the computing grid as shown in Figure 53 is accomplished 
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Figure 51, Rectilinear grid representing the annulus; 
location of hydrographs shown by (0). 
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Figure 52. Polar grid representing the annulus; 
location of hydrographs shown by (0). 
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Computing grid for polar system representation of the annulus. 


Figure 53. 
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where AS* = 21.05 kilometers and AT* = 19.65 kilometers. Approxi- 
mately the same number of computational points are employed in this 
system as that of Figure 51. The actual area of the annulus is well 
represented in both computing grids. Although the inner and outer 
boundaries are straight-line segments in the numerical grid sense, 
the portrayal of the annulus is certainly more representative than 
that shown in Figure 51. The positions of nine hydrographs in 

Figure 52 (denoted by large dots) correspond to the hydrograph posi- 
tions in the rectilinear grid system. For the polar system, the 
maximum numerical time step is approximately 600 seconds. As a 
practical note, the stretched shelf coordinate system, although re- 
quiring less sophisticated programing and allowing better physical 
and mathematical portrayal ot the boundaries, generally results in 

a smaller numerical time step than a rectilinear system of comparable 
scale. 


2. Results. 


For comparison purposes, a time step of 540 seconds is used in 
both grid systems for the numerical solution. As noted above, the 
polar grid program should compile and execute faster on the computer 
than the rectilinear grid program. On the IBM 360/65, the polar 
grid program required 0.47 minute to compile and 0.39 minute per 
100 iterations. For the rectilinear system, the program consumed 
0.59 minute to compile and 0.61 minute per 100 iterations. 


The initial surface elevation condition which is imposed on both 
grid systems is that d0H/d0 is a constant. For the rectilinear grid 
system, 


aes 4 =1 fe il 
H@,37,0) = 1 - | tan G3) é (92) 


where only the i,j indices appropriate for the annulus are employed 


in the above relation. For the polar grid system, the same initial 
surface condition is determined as: 
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The initial condition of constant azimuthal slope and no radial vari- 
ation of H along a line of constant 9 implies that there are 
several modes of seiching present; however, the dominate mode is for 
m= 0 andn-= 2. This is indicated by Figures 54 and 55 which show 
the simulated surface topography at 27 hours in the rectilinear and 
polar grid systems, respectively. The solution in the polar grid 

is more representative of the analytical solution than that shown 

in the Cartesian system. 


The hydrograph at 6 = 0° and r = r, is shown in Figure 56 (a) 
where the water elevation as determined from the polar grid is the 
solid line and that from the rectilinear mesh is the dashline. 
Figure 56 (b,c) shows the hydrographs along 6 = 0° forr=r 
(the average radius of the annulus) and r = ro, respectively. The 
three hydrographs along 6 = m/2 and the same r positions are shown 
in Figure 57 (a,b,c) and along the nodal line by Figure 58 (a,b,c). 
The average period of oscillation as determined from the hydrographs 
of Figures 56 and 57 is approximately 26 hours in the polar system 
and 28 hours in the rectilinear grid. The error in the period of 
oscillation (about 8 percent) for the rectilinear system is most 
evident in the figures by noting the lag of the dashline with 
respect to the solid line. The longer period of oscillation is 
directly related to the stairstep boundary. Effectively, the length 
of the basin is increased by the reflections introduced by these 
boundaries. This distortion is more than academic since many recti- 
linear grid models of enclosed irregular bays require adjustments 
to reproduce the fundamental seiching mode. 


The analytical solution at any point is a smooth function of 
time. The solid lines in the various hydrographs portray this 
feature better than the dashlines which are contaminated by high- 
frequency spurious oscillations. The nodal line in the polar grid 
solution as evidenced by Figure 58 (a,b,c) remains fixed at 6 = 1/4 
which agrees with the theory. This is not observed in the recti- 
linear system. Although the hydrographs in the rectilinear grid 
are not positioned exactly on 6 = 1/4 (actually, about 44°), the 
Magnitude of the oscillation about 1/4 is approximately 2 to 4 
times larger than expected. 


SiG 


Figure 54, Computed water surface topography in 0.0l-meter 
contour increments at 27 hours in the rectilinear 
grid system. 
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Figure 55. Computed water surface topography in 0.01-meter 
contour increments at 27 hours in the polar grid 


system. 
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Figure 56. Hydrographs along ®@ = 0° and vy = rq (a), 
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Figure 58. Hydrographs along 6 = 1/4 and r = 19) 
(a), Ym (b) and rp (c) for the polar 
grid (solid) and rectilinear grid 
(dashed). 
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V. STORM SURGE SIMULATION 
i, Jnitrerclceine (Cereller. 


Hurricane Carla was an immense, slow-moving, and meandering storm 
which struck the gulf coast south of Port Lavaca, Texas, at 2000 G.m.t. 
on 11 September 1961. The radius to maximum winds was 50 kilometers. 
The atmospheric pressure drop across the storm was 75 miltlibars. The 
maximum sustained winds were approximately 51 meters per second. 


For numerical stability a time step of 180 seconds was used in the 
computations. The surge simulation was performed for a 66-hour period 
from 1800 G.m.t. on September to 1200 G.m.t on 12 September. The 
boundary conditions are those given in Section III. The wind fields 
are analytically deformed in a manner previously discussed. The ob- 
served (squares) and computed (solid line) water levels at selected 
grid points along the coast are shown in Figure 59. The observed 
water levels have been corrected for the astronomical tide. At the 
start of the computations, the water surface along the coast from 
Sabine Pass to the Matagorda Bay area was elevated by approximately 
1 meter. Consequently, for comparison the computed water levels 
have been raised to correspond to the local sea level condition. The 
water surface topography in 0.5-meter contour increments at 46, 50, 
and 54 hours are shown in Figure 60. The computed maximum coastal 
surge profile corrected for the astronomical tide is presented in 
Figure 61. Indicated on the figure are the observed high-crest values 
from tide gages and water marks. The depression of the water surface 
to the left of the storm track (as viewed from the sea) and the 
low-peak surge in the Aransas Pass area results from the offshore 
wind and from simulating the coast as a wall. It should be emphasized 
that the surge model does not consider the coastal flooding from the 
shelf surge nor includes the influence of the increased communication 
between the semienclosed bays and the open sea. While the computations 
agree fairly well with those observed for the peak coastal surges, 
especially for the area to the right of storm landfall, and for the 
Galveston hydrography, there is some discrepancy for regions far from 
the storm track such as Sabine Pass. This might be attributed to the 
problem of proper deformation of the wind field (both speed and 
direction) to reflect the influence of land. 


Considerable erosion of beaches and adjacent offshore areas along 
the Texas coast occurred during Hurricane Carla (U.S. Army Engineer 
District, Galveston 1962). Figure 62 shows the water velocity 
V = Qgxi + Qr*j)/D in centimeters per second at selected grid points 
for the same times as that in Figure 60. Of practical importance in 
assessing the erosion potential and transport of the material is the 
alongshore current, Qcox/D - In shallow water, the current is uniform 
from the surface to the sea bed. Values of the current in centimeters 
per second as a function of time at selected grid points are presented 
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Figure 59. Observed (squares) and computed (solid) 


water levels for Hurricane Carla. 
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Figure 60. Computed water surface topography in 0.5-meter 
contour increments for Hurricane Carla at 46 (a), 
50 (b), and 54 hours (c). 
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Figure 62. Computed water velocity in centimeters per second 
at selected grid points for Hurricane Carla at 
46 (a), 50 (b), and 54 hours (c). Vector scale 
in centimeters per second. 
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in Figure 63. Also given is total transport in 10°© cubic meters per 
second as a function of time across the transect lines indicated in 
the figure. 


2. Hurricane Camille. 


Hurricane Camille traversed the extensive low-lying marsh area 
of the Mississippi Delta, moving over the shallow waters of Breton 
and Chandeleur Sounds for 4 hours prior to landfall near Bay St. 
Louis, Mississippi, at approximately 0400 G.m.t., 18 August 1969. The 
radius to maximum winds was about 37 kilometers. The atmospheric 
pressure drop across the storm was 90 millibars. The maximum sus- 
tained winds were approximately 60 meters per second with the highest 
wind gusts estimated at 90 meters per second. 


The surge simulation extended over a 48-hour period beginning at 
0000 G.m.t. on 17 August. A time step of 60 seconds was selected for 
numerical stability. The boundary conditions are the same as those 
used in the Carla simulation. The Camille winds are analytically 
deformed as shown in Figures 43 through 45. The bottom drag coeffi- 
cient was increased to 5 x 1073 at computational grid points in the 
area bounded by 6 < i < 11 and 14 < j < JM to simulate the larger 
resistance associated with the low flooded marsh lying northeast of 
the simulated Mississippi Delta (see Figures 28 and 65). 


Figure 64 shows the observed and computed water levels at se- 
lected grid points along the coast. At Pascagoula, Mississippi, the 
water level increased from 2.2 to 4.5 meters in a 2-hour period 
indicating a strong convergence. In fact, the simulated hydrography 
at the approximate location of the maximum surge (Bay St. Louis) 
showed a similar convergence with the sea level rising from 2.9 to 
7.0 meters in the same time period (figure not presented). The water 
surface topography in l-meter contour increments at selected times 
from 24 to 36 hours illustrating the surge development is shown in 
Figures 65 through 73. The effect of the delta, modeled as a wall 
protruding from the coast, is visible in these figures by the presence 
of high water along the eastern side of the wall. The track of the 
storm with respect to the delta is such that the circulation is 
inhibited in this area resulting in the setup. The water velocity, 

V , at the selected times and grid points as shown in Figure 74 are 
in accord with these observations. The computed maximum coastal 

surge envelope east of the delta presented in Figure 75 shows that 
the calculated peak surges are approximately that observed in the 
Gulfport to Pascagoula area (U.S. Army Engineer District, New Orleans, 
1970). A test run employing a constant wind-stress coefficient, K , 
of 3 x 10°© with the same boundary conditions produced a similar 
maximum surge envelope. There is qualitative agreement between the 
observed and predicted surge envelopes. This is shown by the rela- 
tively low water elevation of 2.7 meters near Hopedale, Louisiana, 
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Observed (squares) and computed (solid) 
water levels for Hurricane Camille. 
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in the marsh area south of Lake Borgne, and a high water mark of 4.8 
meters along the eastern side of the Mississippi River levee. A 
part of the discrepancy between the observed and computed surges re- 
sults from not considering the flooding of low coastal areas and the 
open communication of Mississippi Sound with Lake Borgne which is, 
in turn, connected with Lake Pontchartrain. The tide gage at Shell 
Beach located on the southern end of Lake Borgne crested at 3.4 
meters. Another test simulation was conducted with the hurricane 
track moved 18 kilometers to the east. The results showed greater 
seaward circulation in Chandeleur Sound reducing the setup along the 
levee by approximately 1.5 meters. The location of the peak surge 
was shifted to the east but the surge magnitude was not reduced. 


3. Hurricane Gracie. 


Hurricane Gracie crossed the coast near Beaufort, South Carolina, 
at 1600 G.m.t. 29 September 1959. The radius to maximum winds was 
25 kilometers. The atmospheric pressure in the eye was 950 millibars. 
The maximum sustained winds were approximately 50 meters per second. 


A Cartesian grid for simulating the surge caused by Gracie is 
shown in Figure 76. The reach and grid spacing of this system are 
comparable to the curvilinear grid shown in Figure 32. The thin 
dashlines are a part of the transform-generated coast and seaward 
boundary curve (180-meter depth contour). The surge computations are 
‘performed with both grid systems using the same boundary conditions 
and forcing functions. The stairstep boundary representing the coast 
is the heavy solid line in Figure 76. Another difference between the 
models other than the coastline representation is that the Cartesian 
system's boundary, j = 1, extends farther seaward into deeper water 
than that of the curvilinear system. 


For numerical stability a time step of 150 seconds was used in 
the computations with the curvilinear grid. The time step in the 
rectilinear system was taken as 75 seconds. The surge simulation was 
performed for a 48-hour period with zero hour corresponding to 1200 
G.m.t. on 28 September. 


The observed and computed water levels along the coast in the 
curvilinear (solid) and rectilinear (dashed) systems are shown in 
Figure 77. The inserts on Figure 76 show the simulated hydrograph 
locations. Figures 78 and 79 present the water surface topography in 
the curvilinear and rectilinear grids, respectively, at the same times 
and for the same contour increments. The depression of the water 
surface to the left of the storm track exhibited in both figures is 
Similar to that from the Carla simulation. At the Savannah River 
entrance, the observed water level corrected for the astronomical 
tide is depressed after 27 hours reflecting the offshore-directed 
wind. The effect on the coastal surges and water surface topographies 
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Figure 77. Observed (squares) and computed water levels 


for Hurricane Gracie in the curvilinear (solid) 
and rectilinear (dashline) grid systems. 
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Figure 78. Computed water surface topography in 0.5-meter con- 
tour increments for Hurricane Gracie at 24 (a), 
28(b), and 32 hours (c) in the curvilinear grid. 
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Figure 79. Computed water surface topography in 0.5-meter con- 
tour increments for Hurricane Gracie at 24 (a), 
28 (b), and 32 hours (c) , in the rectilinear grid. 
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resulting from the stairstep portrayal of the coast boundary are 
shown through the comparison of Figures 78 and 79. The curvilinear 
model shows the maximum coastal water level at the time of landfall 
of approximately 3 meters while the Cartesian system predicts 4 
meters. Furthermore, the peak surge of 4.5 meters in the Cartesian 
system occurred at the corner point (22, 19) over 2 hours after 
storm landfall. 


The water velocity, V , from the curvilinear system is presented 
in Figure 80. The figure also shows the development of a cyclonic 
current in the region from Cape Romain to Cape Fear resulting from 
the curved coastline and the extensive Frying Pan Shoals off Cape 
Fear. A similar current pattern is not observed in the rectilinear 
system. 


The computed high water envelopes from both grids are shown in 
Figure 81. The stairstep rectilinear coastline is projected upon the 
curvilinear coastline providing a common basis for the comparison. 
Both systems predict a maximum surge comparable to that of the near- 
est observed high water mark. The effect of the rectilinear coast- 
line on the surge envelope is seen in the figure by the relative 
water level maximums at the interior corner grid points. The pri- 
mary contrast of the coastal envelope of computed surges in the two 
different grid systems is the greater variability of that from the 
rectangular grid. This leads to spurious results as can be seen 
particular for Savannah River. 
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Figure 80. Computed water velocity in centimeters per second at 
: selected grid points for Hurricane Gracie at 24 (a), 
28 (b), and 32 hours (c) in the curvilinear grid. 
See Figure 5 for identification of capes referred to 
in the text. 
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VI. SUMMARY AND CONCLUSIONS 


The development of the conformal mapping problem was completed 
for transforming the interior region bounded by two arbitrarily 
shaped curves and two parallel lateral boundaries into a rectangle 
in the image plane. Specifically, the coastline and seaward boundary 
curve taken as the 180-meter depth contour are transformed into the 
image plane as constant values of n . The bicurve fitting procedure 
was successively tested on five regions of the Continental Shelf in 
the Gulf of Mexico and east coast of the United States. The distinct 
advantages of the orthogonal curvilinear coordinate system over 
rectilinear representations of the coast are the accuracy and ease 
of specification of the boundary conditions. For long wave calcula- 
tions, the numerical simulation of a free long wavelength gravity 
wave in a quarter annulus demonstrated the superiority of modeling 
the boundaries in curvilinear coordinates. 


A two-dimensional time-dependent numerical storm surge model 
employing orthogonal curvilinear coordinates is presented, Since 
the transformation is conformal, the associated modification of the 
equations is minimized. The final coordinates allow for the greatest 
resolution near the coast over the area of principal storm surge 
development and modification. A second transformation is performed 
independently on each of the curvilinear coordinates to provide the 
desired spatial resolution in the shallower region near hurricane 
landfall. 


Results of the numerical simulation in transformed coordinates 
of the storm surges along the gulf coast induced by Hurricanes Carla, 
Camille, and Gracie which crossed the east coast are presented. The 
surface wind fields are analytically represented. Additionally, a 
systematic procedure is employed to deform the symmetric wind fields 
of Carla and Camille in the nearshore region, The deformed Carla 
wind representation produces better agreement between the computed 
surges and that observed than a similar simulation with a symmetric 
wind. Although the pattern of windspeeds produced by the deformation 
procedure is in good agreement with those provided by NWS, the 
problem of specifying the proper orientation of the wind-stress 
vector in the nearshore region requires further study. 


The surge simulations indicate that the model produces results 
in good agreement with the observed peak surges and hydrographs, 
especially in the area to the right of landfall. Moreover, the 
results from the simulation of the Hurricane Gracie surge indicate 
superior rendition in the curvilinear grid compared with that of the 
Cartesian grid. Certainly one contributor to any discrepancy between 
computed and observed water levels may result from the portrayal of 
the coast as a wall. Proper inclusion of attached bays, lakes, and 
flooding of low-lying areas should be pursued, 
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APPENDIX B 


CONFORMAL MAPPING COEFFICIENTS 
FOR THE 


EAST COAST REGION FROM CAPE KENNEDY TO PAMLICO SOUND 
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Note: All quantities are in x,y units. 
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—06.541 89CBE-O1 
OeSSSE3SIOE-Ol 
—-Ocelti3sss5ie-o1 
—06«3498875E-02 
Oc2114é711F5-01 
—-Oe1866217E-O01 
—-02e5141602E—-02 
Ool?i24i4E-O1 
—00713S5S786E-02 
0e3350719E-O02 
Ov 1991S504E—-02 
—004495025E-02 
Oc22432E5E- C2 
0e2640530E-02 
-062747S52E-02 
Oc 361 95S28E-02 
—029984830E-03 
— Ce 987S367E-03 
0«e2608439E-03 
06.1251 844E-02 
02404 8S0B8E-03 
—0 e858 336SE—03 
06326084 9E—-03 
Oe867TSSS6E-02 
—O0e94698S1E-03 
OeS6455SSE- 03 
O0e2324108E-03 
02161 OS@3E-03 
—0¢6130417F-03 
0e5053631E-03 
Oo278SG015E—-03 
—0 04364 226E-03 
-0e3461981E-05 
0e394108SE—63 
—002265753E-03 
Oe 637 72 C4E— 04 
0 e397S5S330E—-04 
0e38307732E—-04 
—-0e106580C6E-C3 
O0e4872639 E-04 
065831S511E-—04 
0640531 38E-04 


100 
101 
102 
103 
104 
105 
106 
107 
108 
ic9 
110 
1i1 
112 
113 
114 


B 
n 
—0e 656451 6E-04 
—0e1748351E-04 
—O«e 3665496F-04 
Oo 10043 76E— 03 
—0 e 6544378E-04 
—O02 S$S63962E-05 
0e2498129E-04 
—0e4729294E-05 
—Qe2e £S413S0E-04 
Oe 330S452E-04 
O. 2570346E- 04 
—02e1S589919E-04 
—-0466349276E-04 
Oo £02411 7E-04 
068936679E-05 
—0.2S8S5104E- 04 
0e1793947E-OS 
0. 28843 07E—-04 
—-0e228911556-04 
—-0e16792i0E-04 
021770000E~-04 
0.2.1685304E-04 
—O20e 267449SE—-04 
0247G2071E-O05S 
06 f7401914E-05 
—O2e5577534E-05 
—022972263E-05 
Oo 7447946E-05 
~0.6316889SE-O5 
-00135627SE-05 
—-O5« 2609007E-OF 
0e367S5S452E-05 
—-02¢1076716E-O€ 
0.£5550138E-06 
—00e 419292 8E- 06 
—00e1196023E-O05 
—-0.834S894E-06 
Oc 311S5964E-O5S 
—O0606 21675S5S6E-05 
—0.e 1715321E-06 
Oe SITO421E- OE 
021282769E-06 
—Oe 213C69SE-O058 
021366156E-05 
Oc 11289S9E—05 
—O02e 1028800E-O5 
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C 
n 
—006221421E-04 
~0.69761238E-05 
—00147282S5E-04 
O.o9884777E— 04 
—0 eS9B7T761E-04 
00511 286SE—06 
063470853E- 04 
—009627895E-05 
~0.0466117SE-04 
0 «3855478E-04 
023692045E-04 
—021510879E—04 
—005784341E-04 
0e508872 0E—04 
001105513E-04 
— 02 294S38SE- C4 
0 e5432613E-05 
002892528E—-04 
—002226605E-C4 
—00¢1523216E£-04 
001914123E-04 
001635740E-04 
—0.2516 768E-04 
0eSO086707E-05 
0e8127184E-05 
= 065464967E-05 
—0«2235201E-05 
007659644 E-05 
— 0¢2746120E- 05 
—061211409£-05 
— 02 2086289E-C5 
023609181 E-05 
— 0« 6883366E- 06 
0e6827090E-06 
—0.1 706489E-06 
—001206698E-05 
—005846961E-06 
0. 3220246E-05 
—0.2048756E- 05 
—001 886759E-06 
0o1148447E-05 
0 013475216—06 
— 0. 2057237E-05 
00140652S5E-CS 
00123400SE-05 
— 0.101 8ES1E- 05 


B 
n 


-001313448E-05 
Oe. 1556963E-O§5 
-0 6 8436349E-08 
—O5c 752393 8E—06 
~001622793E- 0€ 
021139069E-05 
—0« 652760 2E-06 
—0 6 7496735E-06 
Oe 7258587E-06 
Oe S221784E- CE 
-0 2.961 2895E-06 
Oe 200901 OF-06 
0 6 3099258E-06 
—06. 1067625E-906 
—0. 2802620E- 06 
00 2269661E-06 
Oo S356683E- 07 
-0 61395750E-06 
—O5. 1529258E—06 
Ce Z215562SE- 06 
-00e 3408843E-07 
—002487944E-07 
-007241448E- 07 
0. 6702049E-07 
Oe 1388121E-07 
—-0e4186880E-07 
—O. 3S01641E-07 
021001300E-06 
-0¢4257138E-07 
—0.4805924E-07 
0e71350S7E-08 
Oo 74476 08E- 07 
-065374578E-07 
-0e2760153E-07 
Oe 3739260E-07 
0 .61614653E-07 
—0. 5449028E-07 
0¢2701048E- 07 
001567296E-07 
—0201215744E-07 
~0 6 2149563E-07 
0. 1988487E-07 
Oe 3930207E-08 
-0.7725944E-08 
— 0. S216542E-08 
0.61S577168E-07 


14] 


C 
n 


-0 061246 324E-05 
0.6 1568S67E—05 
0 0667168S57E-07 

—007546471 E-06 

—~001113204E-C6 
0 0115C6875E-05 

~0. 616 8775E—-06 

—0e75°270S5E- 06 
067S5807S50E-06 
005233 002E- 06 

—0 0e9294666E-06 
00 2025572E-06 
0¢3316641E-C6 

—001017822E-06 

— 00 2645621E—06 
0 0230 74S0E-06 
001065546E-06 

—001389416E-C& 

—001394174E-06 
0021681 S3E— 06 

—0 e2498253E-07 

—002277810E-07 

—0.6482588E- 07 
006857 S88E-07 
0c 2008024E—07 

-0e46101790E-07 

—003311435E£-07 
0010054S8E- 06 

—0-3833354E-07 

— 00471 6878E-07 
021083663E-07 
007493 753E-07 

—0e5108661E-07 

-022660881E-07 
Oe 39585S93E-07 
061627980E— 07 

-0252043995-07 
0027223 71E-07 
Oel 71 SO007E—07 

—0.1168 664E-—OF 

—00201 88EEE- 07 

'Qc2019971E-07 
0. 4938¢83E—08 

—0067394348E- 08 

—0«68195819E-08 
001577923E-07 


EAST COAST REGION FROM PAMLICO SOUND TO PENOBSCOT BAY 


Sie 


OMNAHWPUD 


Note: 


APPENDIX C 


CONFORMAL MAPPING COEFFICIENTS 


FOR THE 


3s 


w 
it 


B 
n 
O«e 8453177E O1 
—-O56e2573375E OZ 
O022456377E O1 
—-0¢4028893E O1 
—02e1294892E O1 
Oe 7517704E OO 
—-05.26833593E OC 
0e8799276E O00 
—0.e 568954 8E 00 
0264709852E 00 
-O50 350S102E 00 
—-00eSS0869S6E-01 
Oell06679E OO 
—O0« 3838463E-01 
—O0« 2045463E-01 
-022993307E-01 
0655191 73E-01 
-—0 e3588944E-01 
OO. 2528971E-O1l 
—O00« 2555390E- 01 
0e2158327£-01 
—06« 2883616E-O1l 
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Oe 21852463E 02 


Ce1l045535E 03 


37420 


C 
n 
—Oe7390710E 02 
Ooll466E56E 02 
—-0e131672S9SE 02 
Oe8089012E O1 
—-0eS803839F O01 
00e«2271722E O1 
—-0e6068285E 00 
0e8849853E 00 
—-00e1988623E O00 
Q023153S9S88E 00 
—-0021405933E 00 
—-020«16112S7E OO 
Oe20470SSE CO 
—-0e7570612E-01 
0e3443124E-02 
—0 e4593338E-01 
Co 9201282E-01 
—-00e688 24S535-01 
00e5296847E-01 
—-00e439624EE-O01 
G0e3201135E-01 
~—0e3150356E-01 


All quantities are in x,y units. 


B 
n 
0c1704723E-01 
—Oe 821 84 825-02 
0e5179290E~ 02 
-02301S571SE-02 
Oc 793892S5E-02 
-0 e 7356405E-02 
O. 21203 O3E-02 
—062383103E- 02 
02.1206459E-02 
—00 4140623E~ 03 
0e3749S5S28E-02 
O5oc 4151922E-03 
—O«. 624251 BE- 03 
Oe45S18004E-0 23 
-02e1179293E-03 
—-0e110SO91E-02 
—O50 261 8267E-03 
Oe 33287 O8E- O02 
—-0 e 219220SE-03 
Oe 1223145E—- 03 
—-028201750E- 04 
Oc 1276652E—-0 3 
—-0. 740803 7E-904 
0.402277T1E-04 
—05. €5351S50E-04 
02€436008E-O5 
—0e7533055E-05 
0. 208 7624E- 04 
—-02 8%48817E-OS 
—O0 273 7144E—-05 
02 109S4GEE-04 
-0.65011S9E-05 
—O50« £406307E-O5 
—-0 2e2243007E-0S 
Oe 4238569E-05 
—-0e1688697E-O05 
0263522237E-05 
—00. 2636244E-905 
0e1468087E-0S 
—-Oe 1664401E—-905 
00 16220326-05 
—-0e1230003E-05 
—0.222695T1E-—06 
0e&6773S8E-07 
0. T2Z02404E-06 
—O0« 679648 3E- O€ 
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C 
n 
0.61959118E-O1 
—00e7724181 E-02 
Oe 3032119E-02 
0 01962418E-03 
Co 488653 0E-62 
—0e41S50771E-02 
021651054 E-03 
—00e11456i6E-02 
0 e4586522E-03 
Ool 666064E-—03 
00132 8797E-C3 
0e599561S5E-03 
—0e6716528E- 03 
0045008S59E- 03 
-0°e8762798E-04 
—0e1055562E-03 
—002268526E-03 
Oe 295065 BE— 03 
-0.16891S6E-03 
0ce8370787E—04 
—02382705SE—- C4 
006963981 3E—-04 
—0c4502C0S0E-04 
0e2078711E—-04 
—0e5105330E-04 
0.1149840E-05 
—00420 369S5E-05 
0. 3067S82E—04 
—0208395030E-05 
—062279463 E-05 
Oo 9629463E-05 
—004204831E-05 
— Ceo 72332S0E-C5 
—0e388391SE- 06 
002996593E-05 
— 06521 7331E-06 
0<«2865336E-05 
—022091638E-C5S 
O0.1271E€15E-O05 
—061547839E-05 
0c 166901 CE-05 
—-0.1246244E-05 
—001S536769E-06 
0e4681078E-—07 
0«e7948224E-06 
— 0. 727 0842E- C6 


1co 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
112 
112 
113 
114 


B 
n 
0064362259E-06 
—O5«. 1076482E-06 
0e24S59019E-07 
—O6 8354027E-07 
—O«. €96S583E-07 
GeS715553E-—-O0S 
—-02.10971125-O07 
Oe 7904828E-07 
O062604376E-07 
—-00 481 7723E-07 
-Oe.lO023679E-07 
O52. 35211 895-07 
—-O0eS5271849E-07 
O22211178E-O7 
—0O. 7624720E-08 
0e1252196E-07 
—0e 4G695075E-—08 
00 €09S036E- 08 
—0e6119723E-08 
Oe €175922E- 0S 
-Oe 3125592E-08 
Oo. 560471 8E-08 
—0e6158700EFE-08 
02e1393683E-08 
Oe 2143876E- 0€ 
-0 e3371285E-09 
—0. 201i 3834E-0S 
Oxo i21930R8E- OS 
-02e4478S535E-09 
—06 26221 87E-O09 
0e1349961F£-09 
—0.8012752E-10 
0. 1434S38E- 0S 
—-0e 7944631£-10 
0e469S801E-09 
—0 9 3936562E-09 
—0.66882636E-10 
O0« Z980914ZE-16 
0.6826086E-10 
—Oe«1151273E-0S9 
OelO01L11LO2SE-0S 
-O56. 3327474E-10 
Oc 243134S5E-1¢€ 
—0e¢8947284E-11 
—05. S061 761E-12 
—0-.1896980E-10 
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C 
n 
0 e4967463E-06 
—0014C8421E-06 
0066843044E—-C7 
—-0c¢109C0609E—-06 
—024082S53E-07 
—-0cei3E0970E-07 
O0e3200026E-08 
Oe /S80621E-07 
0230SS97T9E-07 
—CoS 74 C57SE-07 
—-0e954938SE—-08 
Oe3772427E—-O07 
—-0654265S9E-07 
0e2568847E-O07 
—0269943246E-08 
001546062E- 07 
—-00e«6723383E-08 
00e6111858E-908 
—O06e725S402E-08 
001 844SC0E—-08 
—-023619169E-08 
02605S5163E—-08 
—-0e6123S51E—-08 
0¢61397110E—-08 
0e2388478E—-08 
—-Oce48S838 7E-C9 
—-0e35220485-10 
- 0.05434 760E-10 
—0 62261 494E-909 
—-0.2368781SE-—09 
0e2747322E-69 
—001470834E-09 
062115 7SEE—0G9 
—-O0e977741S9SE—-10 
0e4927398E-09 
—-0.23841 8ESE-09 
—00e7195172E-10 
Oe47375235E-10 
0e5399223E-10 
—-0069696301E-10 
0e-9G092571E-10 
—0c«196S5071E-10 
Ooi71162SE—-10 
—068525016E-12 
—023881352E-11 
—061490SSSE-10 


B 
n 
—O2 7091S71E-11 
Oe 1821453E-10 
—-0¢ 3987485E-11 
—02 $260623E~-11 
Q5e7307409E-11 
Oc 94B87284E-11 
—Oe 1281272E-1C 
O02eS417571E-11 
—05. 4854271E-11 
O021272S995E-11 
—0.8162463E-12 
0. 2438756E-11 
—0201364847E-11 
Oo. 7146806E-12 
—-0e4318938E-12 
GO.1000S46E-11 
—0.« 2052629F-11 
0e615S0933E-12 
Oo 4740078E-12 
—O5« 39357S54E-1ée2 
0061734885E-12 
Oo 122661 3E-12 
—-021624964E-12 
0. 237601 8E-13 
02¢1584382E-12 
-007575429E-13 
—0.160i1686E-13 
-0262932933E-13 
Oo 14569 78E-12 
—-O0211467427E-12 
062950459E-13 
Oe 2668 76i1E-13 
—0.2074897E-12 
—Oe 2806369E-13 
002E75362E-12 
00247861 3E-13 
0. 1650332E-13 
—-021226408E-14 
Oo TO71394E-14 
—0O.S688173E-14 
—0¢ 3497994E-15 
Oe 4112430E-14 
—-0e3376250FE-14 
—0.3620927E-14 
0651592S5S7E-14 
-0264390035-15S 
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C 
n 
—O0e7386E324E-11 
0061942668E-10 
-O0e311S7ESE-11 
=—Oe 95221 16E—1 1 
OeS7TO1L047E-11 
O0.874S027E-11 
—0201152168E-10 
0047035i9E-11 
—00e587 0245E-11 
O0e5C34556E-12 
—002376758E-12 
Ce2222751E-11 
—-02e1078S58E-11 
Oe6671415E-12 
“On S27 SMyyse = Ue2 
0e1055302E-11 
—-0e2062S67E~-11 
Oe?%127636E-12 
0<428839E9E-12 
Wo SOS GOSS ie 
06i294625E—-12 
OcolBSSBOLE-12 
—-00l192G000E-i12 
0e6448195E-13 
063484 010E-14 
—-0.5707148E-13 
—-Oelé698161E-13 
—002245293E~-13 
Oeil 4S98009E-12 
—0e1142878E-12 
0e3582161E-13 
0e2435500E-13 
—-0614775VY2E-13 
—-00e30550S6E-13 
0.31 08SE5E—-13 
—022643752E-13 
Ocl9S33177E-13 
—~Oe2t20201E-14 
Oe 864977S5SE-14 
—-Co 9941 355E-14 
Oe3751107E-15 
Oe4258614E-14 
—-Oe3143469E-14 
—063506923E-14 
0.5141 81SE-14 
—-Oe3103656E-15 


B 
n 
-O4. 1SS0099E-14 
021889396E-14 
—-0e1462261E-14 
Oe SSSIOISE-16 
-O021863618E-15 
Oo 4924783E-15 
—Ce€9751L17E-15 
022073116E-15 
OG. 2038455F-15 
022677031E-15 
-026432034E-15 
00e4125446F- 15 
-021303287E-15 
—-021122934E-15 
0-e4706S5S50E-16 
Oo 55621126-16 
—-Oo 728751 7E-1€ 
0e6455897E-16 
—-0e1810812E-17 
—-0e25194S57E-16 
—O5« 2583353F-16 
0eS397427E-17 
0e284996S5E-16 
— Oc 435331 7E-16 
00e2342945E-16 
Oe 6163340E-17 
-0.21038789E-16 
—02254589S5E-17 
Oe 2582216E-17 
—-0212928745-16 
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C 
n 
—O0el1644377E-14 
0621 769S3E- 14 
-0 e15E2678E-14 
Oc 3085169E~-15 
SWS Zoos SSVS 1S) 
0e5232375E-15 
—0067213222E-15 
Oe2703232E-15 
Oe2139917E-15 
Oe2874004E-15 
-O02e«6162772E-15 
004083402E-15 
—-069925271E-16 
—-O0«0l216983E-15 
0e7034376E-16 
024932992E-16 
—0e5844C15E-16 
0e6172391E-16 
Ooe71 76773E-17 
—002643164E-16 
—0e 2060502E-16 
06. 9356655E-17 
0 03153878E- 16 
—-00e4277073E-16 
0.2463311E-16 
O0«7743982E-17 
—-021057159E-16 
—-0e5732719F-18 
0e8081919E-17 
—-0e11264SEE-15 


GULF COAST REGION FROM LAGUNA MADRE TO MARSH ISLAND 


=) 


ODNANE UA = 


Note: 


APPENDIX D 


CONFORMAL MAPPING COEFFICIENTS 


FOR THE 


0262110953E 02 


B = Ce9SO7TR64LE O02 


= 360.0 


B 
nh 
—004349387E 02 
—O50o Z2166200E OZ 
0.1242251E 02 
—-0.5063650E 01 
0e1671185E O01 
—-O« 3597688E O1 
0.1419286E O1 
—Ocol1179494E 01 
0.4119087E OO 
—O003775935E OC 
0e7971597E-0O1 
—0. 1080030E 00 
—0e1290084E-01 
—00 SS63697E-01 
—O«¢ 1379531E-01 
—0.3508482E-01 
Oe 6604981E-02 
—001430117E-01 
06 93927S5S7E-02 
—0.1132046E- C1 
—001077687E-02 
—02¢1274931E-01 


All quantities are in 
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Cc 

n 
—-Oe6817682E O02 
001942726E 02 
—-Oell73821E O1 


~=-001447072E-01 


—-001921633E O1 
—066733062E OO 
OcelS66718E C1 
—-06309135SE 00 
Ocl335779E O00 
004832629E~-02 
0 eS321220E-01 
Oe1E0S034E 00 
-00215S98e8e8E 00 
Oe14032SSE 00 
— 0c 6387305E- 01 
0 096i 8145E-01 
—004127430E-01 
0¢213399S5E-O1 
—0.2286828E-O01 
Qce23924E0E-01 
—0 e2269067E-02 
06632 0797E-—O02 


x,y units. 


68 


B 
n 
Oe 4066765E-02 
—0O.2 L667407E-02 
0e3995597E-0 3 
~Qe 188621 9E-02 
0.2907 7224E-03 
—Ooe 155919 8E-02 
0eGS190331E- 02 
—-02114382S5E-02 
O56 37S5230E-03 
—-0e6712531E-02 
Oo 4371286E-03 
-0. 583528 e8E- 02 
02.1835159E-03 
—06. 251 2983E-03 
—-O0210435S64E-04 
Oe 1022287E-04 
—0-¢ E843900E-O0S 
—O00e127881 3E-04 
—050e 261 3426E- 04 
—-02¢1440295E-04 
—-0e5079823E-05 
—-O0e 238344 7E- 04 
022685109E-04 
—-06.4311551E-04 
0 «26351 92E-04 
—06 2369922E-04 
Oe 1060475E- 04 
—02e11909G66E-04 
O-. #96911 9E-05 
—-029435032E-05 
Oo 3592144E-905 
—0-3163786E- OF 
0.8069733E-06 


=O. 31561 97E-05 
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APPENDIX F 


NUMERICAL ANALOGS OF SURGE EQUATIONS 


Admissible finite difference analogs of equations (40) and (41) 
which are nominally centered at i, j , n (bearing in mind the 
storage of variables in Figure 36), are as follows: 


[Qa (isd mt1) - Qo, (G,j,m-1]/2at 

- £[Qpe (4.5.41) + Qpe(i,j.n-1)]/2 

+ (gD/F(i,j)u(i)) [HG+1,j,n) - Hy (i+1,j,n) 
= GF on) < Gol jaja) |//ans 


= to, (isj,n) - (K,Q(i,j,n-1)/D7)Q,, (i,j n+) (F-1) 
and 


[Qpe(i.5.n*1) - Qp,(i,i.n-1)]/2at 
Oe Gi and) Oey Gai n=t)))/2 

+ (gD/F(i,j)v(j)) HG, j+1,n) - Hy (i,j+1,n) 
- H(i,j-1l,n) + H, (i, j-1,n)]/2aT* 


Some) COG jyn-1)/D Ore Ej net) (F-2) 


where Q is as defined by equation (54) and D as defined by 
equation (61) is the arithmetic average of the four values of D 
about the point i,j at which the flow is evaluated. All terms are 
spatially centered at i,j and all but the bottom friction terms 
are exactly centered at time level n . The latter involve flow com- 
ponents at the new level n+l and at the old level n-1 within the 
flow magnitude term Q . This form is known to lead to a stable 
algorithm (Reid and Bodine, 1968). The Coriolis terms are exactly 
centered at n as well as i,j . Equations (F-1) and (F-2) repre- 
sent two equations in the two new Q components Qcox(i,j,ntl) and 
Qrx(4,j,n+1) in terms of quantities at previous times and are 
readily solved for each component individually. The resulting 
explicit relations for these components are given by equation (55) 
where G; , Go , G3 are defined in equations (58) to (60). 


The centered finite difference analog of the continuity equa- 


tion (42) leads directly to. the explicit relation for H at the new 
time level given by equation (57). 
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APPENDIX G 
MODEL VERIFICATION 


Prior to the surge computations the conservative properties of 
the numerical model were examined. The model was applied first to 
seiches in an enclosed square basin 12 kilometers on a side and 5 
meters deep. The number of grid increments along each side was 20. 
To invoke the most severe seiching mode, the water surface was 
initially assigned a linear slope downward from one end of the basin 
(j=JM) to the opposite end (j=l). The initial setup for the nonro- 
tating basin computations was 50 centimeters, but this was reduced 
to 10 centimeters when rotation was allowed to suppress unnecessary 
fluctuations. The time step taken for all test runs was 50 seconds. 


The simulation of a uninodal seiche in a nonrotating basin was 
considered first. With the specified time step, tne period of the 
severest mode is theoretically 68.57 At. 


Figure G-1 shows that the computed fundamental period is within 
one-time increment of the theoretical period. 


Water Elevation (m) 


! 
0 30 60 90 120 150 180 210 


Time (nAt) 


Figure G-1. Computed water level at grid point (21,20) 
for uninodal seiche test case. 


The time variation of the volume, potential, kinetic, and total 
energies are given in Figure G-2. Each of these quantities are area- 
weighted sums. This evaluation scheme should exhibit small perturba- 
tions more readily than unweighted sums of the variables over the grid. 
The variation of the volume, potential, kinetic, and total energies is 
small, with the potential energy showing the largest deviation. 
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Energy in relative units 


Potential energy (P) 
Kinetic energy (K) 


Figure G-2. 


energy (T) 
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Volume, potential, kinetic, and total energies for 
uninodal seich test case. 


157 


210 


Rotational effects on the seiching mode were checked by specify- 
ing the Coriolis parameter as 0.75 05, 1.00 og, and 1.25 o9 where oo 
is the angular: frequency of the fundamental mode. The computed fre- 
quencies, o , were compared with those obtained by Rao (1966) and 
Platzman (1972). 


Table G-1 gives the comparison with Rao's and Platzman's results. 


Table G-1. Comparison of computed values of c/o, with those obtained 
by Rao (1966) and Platzman (1972). 


5/59 
SS S=SSSS—EEE Average 
From Rao From Platzman computed 0 
£/o9 (1966) (1972) 0/09 oI ee 
0.75 0.769 0.76 ed 
1.00 ORZS OQ. 72il ORS 4.0 
WEBS 0.686 0.683 0.71 3.9 


Although the errors occurring with f/o, > 1.00 are large, these 
values of f/o, represent stern tests of the algorithm. It is 
definitely encouraging that no evidence of instability is indicated 


‘in the simulated hydrographs. 
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APPENDIX H 


WIND DEFORMATION PROCEDURE 


Assuming that the wind-stress components (co OD) 5 eo) 
are known in the stretched shelf coordinate system from the applica- 
tion of a symmetric wind field model, the problem in part is to 
identify a region for a given hurricane and coastline where it is 
appropriate to alter these stresses to reflect the influence of land. 
The proposed deformation equations are applied only at points which 
are located within this region. Existing charts of hurricane winds 
and the investigation by Graham and Nunn (1959) provide the basis for 
the empirical deformation formulas. In this manner, the analytical 
representation obviates the detailed input of a massive sequence of 
digitized wind field data that conforms with the observations near 
the coast as well as offshore. 


The wind-stress components, rent) and on oe are 
altered at a point for which Y; > Yp according to the relationship: 


(SYM), 
Tox BF. Tox I, 2 (H-1) 
and 
SYM m 
a ae (H-2) 
where m is a constant chosen to be 2 and I. is an influence 
factor given by: 
= - - R 2 = 
I 1 Del (Yy Yp)/R I 5 (H-3) 


where Rp is a constant taken as the average radius to maximum winds 
of the hurricane, Dg is a distortion factor, Yp is the shortest 
distance the point is from land, and Yj ais the distance the 
influence region extends from the coast relative to the point in 
question (Figure H-1). In general, Yy should depend on Xp , the 
distance the point in question is from a centrally located point along 
the coast (X,,y,.).- This latter point is generally assigned the 
coordinates of intersection of the hurricane track with the coastline. 
From Figure H-1, the distances Xp and Yp are given by: 


1 
Kye Bl @& = a) ® Gas Ye alll ; (H-4) 
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Yy =x, - x)? +, - v2 (H-5) 


where (X5,y,) 1S a point along the coast determined by the inter- 
section of a line normal to the coast and passing through the point 
in question (X,,y,). The sign convention is that Xp is positive 
if the point in question is to the right of the hurricane track and 
negative if it lies on the left. In general, the influence region 
is defined by: 

Yr =[ Fi (Xp) + Sp] Re ; (H-6) 
where Fy; is the influence function of Xp, S_¢ is a shift factor 
which permits the influence line to be moved normal to the coast, 
and Re is a range factor. The influence function attempts to re- 
late the degree of distortion in the wind isovels relative to the 
tract of the hurricane and the shape of the coastline. The extent 
of the land influence is given by Yj , which depends on S¢ and 
Re . These coefficients as well as Dg are considered to be 
related to the distance of the hurricane from the point where it 
crosses the coast. 


Additional investigation is required to establish a generalized 
procedure based on physical principles for deforming the winds in a 
manner which faithfully reproduces the observations. Consequently, 
each storm is parameterized individually to accomplish the task. The 
relationships used for deforming the winds of Hurricane Carla are: 


Fr 5 | X Bs RX, |/7 > (H-7) 
Dee sea (H-8) 
S. =f s5- = tan! (0, /2R)) R, , (H-9) 
1.5 peices D, < -5R,, 
Re = ie is (H-10) 
Be 7 025 (D, /R, +5)?, if D,, > SR, 


where RX. is taken as approximately 37 kilometers and 
1 
fh is 2 = 272 RS 
De meal(Ge sa NAC en ele: (H-11) 


The sign convention is that is positive if the hurricane is over 
water and negative if the storm is over land. For Hurricane Camille, 
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these relations are: 
(X46, )/1-5.5 2£ 5 < -3R, 
BS [2S 5 [/ed 2 He =§i, << (H-12) 


| x Wag » HEI SO 


De = 052 5 (H-13) 

S¢ = 0.5, (H-14) 
and 

Re = i360 (H-15) 


The Hurricane Gracie wind fields were not deformed because the 
winds are presented in digital form. Any deformation obtained would 
be dependent upon a subjective analysis which requires additional 
developmental work. Moreover, the selection of Gracie was made on 
the knowledge that the results would possibly be compared with those 
obtained by other investigators using symmetric wind fields. 
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APPENDIX I 
SYMBOLS AND DEFINITIONS 


Transformation coefficients 
Nondimensional wind-drag coefficient 


Depth of water 


Distortion factor for deforming the wind 
Distance of hurricane from the coast 


Maximum expected depth of water at any point 
within computing grid 


Depth of water relative to mean sea level 


Error function 
Coriolis parameter 


Scale factor associated with the orthogonal 
curvilinear coordinate system 


Influence function for deforming the wind 


Acceleration due to gravity 

Water level relative to mean sea level 
Hydrostatic elevation of the sea surface corres- 
ponding to the departure of the atmospheric 


pressure from a constant 


Influence factor for deforming the wind 


Number of computational grid points along the S* and 
ieecaxXes 


Wave number = 1/A 


Wave number of free wave in the annulus 
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X,Y 


Nondimensional bottom-drag coefficient 


Number of Fourier-type transformation coefficients 


in the series 
Atmospheric pressure 


Far field pressure 
Central pressure of the hurricane 


Volume transport per unit width in the S* and 
T* directions 


Polar coordinates 


Distance from the hurricane center to any point 


Range factor for deforming the wind 


Distance from the hurricane center to the region of 


maximum winds 


Shift factor for deforming the wind 


Distance normal to the seaward boundary and along 
the coast 


Coordinates of the stretched shelf coordinate 
system 


Long wave travel time 

Maximum winds of the hurricane 

Weighting factors for the coast and seaward 
boundary curves 


Windspeed at an elevation of 10 meters above the 
water surface 


Rectilinear coordinates 
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At 


Usv 


Coordinates of landfall of the hurricane 
Position of hurricane center 
A point in the influence region 


A point along the coast determined by the inter- 
section of a line normal to the coast and pass- 
ing through (x, >¥,) 


Coordinates of the coast and seaward boundary 
curves to be mapped 


Extent of the land influence for deforming the 
wind 


Plane of (x,y) 

Extent of n (+) in the c-plane 
Arc length 

Surge algorithm time increment 


The alongshore curvilinear coordinate of the shelf 
coordinate system 


Prorating factors for the length of the transform- 
generated curves 


Half the horizontal extent of the region to be mapped 
in the Z or f-plane 


Functions transforming & and n to the stretched 
shelf coordinate system (S*,T*) 


Ingress angle 


Density of air and water 


Total and nth contribution to the "curvilinearity" 
variance 
Bottom resistance stress divided by in nS Se 


p 
and T* directions W 
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Tox? Tpx 


Wind stress divided by Pw in the S* and T* 
directions 


Angle between the € and x axes 


The offshore curvilinear coordinate of the shelf 
coordinate system (&,n) 


Plane of (&,n) 
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